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Abstract 

Hadronic yields and yield ratios observed in Pb+Pb collisions at the 
SPS energy of 158 GeV per nucleon are known to resemble a thermal 
| equilibrium population at T = 180 ± 10 MeV, also observed in elemen- 

tary e + + e~ to hadron data at LEP. We argue that this is the universal 
consequence of the QCD parton to hadron phase transition populating 
the maximum entropy state. This state is shown to survive the hadronic 
rescattering and expansion phase, freezing in right after hadronization 
. due to the very rapid longitudinal and transverse expansion that is in- 

T-j- | ferred from Bose-Einstein pion correlation analysis of central Pb+Pb 

CN ' collisions. 

in 

\ Recent Lattice-QCD theory predicts the disappearance of the hadronic 

phase of matter once the energy density exceeds a critical value of about 1 to 
1.5 GeV per fm 3 [l], giving rise to a continuous, partially deconfined QCD 
state that is governed by the elementary interaction of quarks and gluons. To 
Qh! recreate this phase in the laboratory one collides heavy nuclei at relativistic 

energy, with the goal of ascertaining the QCD predictions, and to pin down the 
decay point from the partonic to the hadronic phase by obtaining estimates for 
the transition temperature and energy density. The CERN SPS Lead ( 208 Pb) 



X 



beam facility provides for a top energy of 158 GeV per projectile nucleon, 
corresponding to a total internal CM-system energy of about 3.5 TeV, to heat 
and compress the primordial reaction volume. In fact calorimetric data show 
that the average transverse energy density exceeds about 2.5 GeV/ fm 3 in that 
volume, in central Pb + Pb collisions. Moreover the study of J/ty production 
H demonstrates a suppression of the yield in such collisions, characteristic of 
the QCD "Debeye" screening mechanism expected in a deconfined partonic 
medium 

If one tentatively takes for granted such indications of a deconfined state 
at top SPS energy one expects, likewise, to receive signals specific to the bulk 
parton to hadron phase transformation bound to occur once the primordial 
high density state expands and cools toward the critical temperature and en- 
ergy density. In the present letter we shall examine the proposal || |[ [7|, |8j 
that such a signal is indeed provided for by the composition of the bulk hadron 
final state, in terms of the yields and yield ratios of the various mesonic and 
baryonic species. To this end we proceed in two steps. Firstly, a synthe- 
sis is attempted between the observation made in microscopic parton cascade 
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and hadronization models |5], |9|, |10| that the nonperturbative mechanism im- 
plied in the parton to hadron "coalescence" transition leads to a population of 
the hadronic species that is dominated by conservation laws and phase space, 
and, on the other hand, the well known conclusion from thermal equilibrium 
models ]7], [11], [12], §] that the hadron yields resemble a hadro-chemical equi- 
librium population. We shall conclude that this apparent equilibrium state is 
not established by dynamical equilibration amongst the hadronic species; it is 
an outcome, and a direct fingerprint of the hadronization transition. In the 
second, main step we address the question of how such a multihadronic popula- 
tion pattern established right at the transition point can survive the extensive 
rescattering that occurs in the course of subsequent hadronic expansion Jl3| 
in central Pb + Pb collisions. This will be answered by confronting the dy- 
namical expansion time scales derived from Bose-Einstein two pion correlation 
analysis of such collisions |TJJ with estimates of the hadro-chemical relaxation 
time scales [|15|, 16], showing that the "explosive" nature of the expansion (17 



prevents major rearrangement of the hadron population. 

Turning first to the connection between the hadronization mechanism and 
the apparent hadro-chemical equilibrium we refer to the partonic cascade 
model of Ellis and Geiger ||. They analyzed LEP data, e + + e~ — > Z° — > 
hadrons at ^fs = 91 GeV in a quark-gluon transport model based on pertur- 
bative QCD. The cascade ends with quarks and gluons recombining to form 
colour singlets with matching flavour combinations in a statistical coalescence 
process which generates heavy "pre-hadrons" that, in turn, decay into the fi- 
nal hadrons and resonances according to their relative phase space weights. 
No hadronic rescattering is considered. After resonance decay the hadronic 
spectra and production multiplicities are compared to the data finding good 
agreement. The authors report a second order sensitivity of the hadron com- 
position to the details of the assumed parton coalescence concluding that, to 
first order, phase space dominance during the hadronization phase overwhelms 
other influences. The hadronic final state is thus the most probable, maximum 
entropy state. 

Becattini has analyzed [7|, [TT| the same LEP data using a hadro-chemical 
equilibrium model based on a canonical partition function. With fit parameters 
temperature T, reaction volume V and a strangeness undersaturation param- 
eter 7 S he finds good overall agreement with the data using T = 161 MeV 
and 7s = 0.67. The data comprise hadron multiplicities for about 25 species 
ranging from tc to Q and spread over 4 orders of magnitude. As there is essen- 
tially no hadronic rescattering occuring in the reaction the apparent equilib- 
rium state can not have originated from inelastic interactions between the 30 
hadrons created per collision event. I. e. there is no chemistry at all between 
the hadrons, they form no stationary state with dynamically maintained pop- 
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ulation ratios. In accordance with the outcome of the Ellis-Geiger model || 
we conclude that the order seen in the hadronic population is born into it by 
the non-perturbative hadronization mechanism. 



Moving on to SPS Pb + Pb collisions at y/s = 17GeV similar data exist Jl 



for hadron multiplicity ranging from pions to Omegas, and similar models are 
successfully applied to these data. Geiger and Srivastava |23| have extended 
the parton transport model || to such low They report that the dominant 
fraction of the hadron yields (near midrapidity) stems indeed from the parton 



to hadron phase transition. They get a reasonable fit to the NA49 data |T 
concerning proton, pion and kaon yields and rapidity distributions, without 
any consideration of secondary hadronic cascading: a baffling result. 

Becattini et al. |I2j employ a generalized hadro-chemical equilibrium model 
based on deducing a temperature parameter of T = 192 ± 19 MeV for the 
NA49 Pb + Pb data set \TB] which extends from pions to the A/ A ratio only. 
Braun-Munzinger et al. || address a data set for central Pb + Pb combining 
all available NA44, NA49 and WA97 information obtaining a satisfactory fit 
at T = 170 ±11 MeV. We note that the two thermal models mentioned here 
differ in several ways, concerning e. g. the treatment of excluded volume and 
specific approach to strangeness saturation. The question of a final, "optimal" 
formulation that avoids an inflation of free parameters is still open [12|, || [19! • 

We ignore this question placing prime emphasis on the conclusion that 
again both the parton cascade hadronization mechanism and the hadro-chemi- 
cal equilibrium ansatz fit the data, the latter employing "temperatures" that 
are slightly higher than those obtained for elementary e + + e~ collision data. 
The large hadronic system created in Pb + Pb through the hadronization phase 
transition can thus not " cool" the hadronic population. Although in this case 
the idea of re-equilibration due to inelastic hadronic rescattering seems appli- 
cable, the above data and results indicate that this is not occuring: the state 
created at hadronization stays essentially frozen-in throughout expansion, and 
its equilibrium features again appear to have nothing at all to do with inelas- 
tic interactions among the 2500 hadrons |18| in the final expansion. Solving 
this puzzle will help to establish the multi-hadronic population ratios as a 
fingerprint of a parton to hadron phase transition occuring in SPS collisions. 

Hadronic inelasticity is rather weak in a hadronic medium at T = 180 MeV: 
the mean CM kinetic energy of a pair of baryons is only 270 MeV and we 
recall that the average pion multiplicity in p + p collisions at this energy |20| 
is only about 0.2. Only every fifth collision is inelastic. Hadronic equilibration 
times should thus be well above 1 fm/c at this temperature, growing to about 
20 fm/c at T = 120 MeV []nj, the temperature of final hadronic decoupling 



from strong interaction [14 



To appreciate these relaxation time scales we need to know about the ex- 
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pansion time scale in central Pb + Pb collisions. Fig.l shows the time profile of 
pion decoupling |H] from the fireball. The Gaussian with mean 77 = 8 fm/c 
and width At = 4 fm/c results from an analysis of Bose-Einstein pion pair 
correlations in central Pb + Pb by NA49 |14| using a dynamical hadron source 



expansion model developed by Heinz [22]. The time scale in this model be 



gins after the primordial source is formed, irrespective of its content at this 
time (partons or hadrons). This time is reached at complete interpenetration, 
t ~ R(Pb)^Q M ~ 1 fm/c. Then the emission clock starts (r = in Fig. 1), 
correctly so as pions can be emitted from the surface of the fireball right away. 
We see that the mean decoupling (freeze-out) time of pions is r = 17 = 8 fm/c 
from this analysis. This includes the average life-time of the partonic phase. 
In order to arrive at an estimate of the mean hadronization time we reproduce 
in Fig. 2 a result of the partonic cascade model by Geiger and Srivastava f23 



that we quoted above to reproduce global features of central Pb + Pb collisions. 
This shows the time dependence of the average energy density in the vicinity 
of midrapidity as exhibited by partons and hadrons (plus resonances), in a 
cylindrical subvolume of 4 fm transverse radius, i. e. in the interior section of 
the source. The partonic era ends at t — 6 fm/c. We may crudely estimate a 
mean hadronization time of 3 fm/c. 

The average time span between hadron formation and freeze-out may thus 
be estimated to be about 6 fm/c (taking proper account of the fact that r = 
in Fig.l corresponds to t — 1 fm/c in Fig. 2). However, during this time 
interval any freshly formed hadron experiences a dramatic falloff concerning 
density and temperature in its local, co-moving environment cell that partic- 
ipates in the overall expansion of the system. A closer look at Fig. 2 already 
reveals the presence of this expansion: the (energy) density of hadrons and 
resonances in the inner subvolume considered here stays constant through- 
out the hadronization phase although it is constantly replentished by further, 
newly created hadrons. This implies that the rate of hadronization inside the 
cylinder must be compensated by the rate of hadrons escaping, thus prevent- 
ing a pile up to unreasonable central density beyond 0.6 GeV / fm? . A more 
quantitative consideration shows that the newly formed hadrons must have 
left this inner subvolume of maximal density and temperature after about 
1.5 fm/c on average, much smaller than the estimated typical hadronic relax- 
ation time of about 4 fm/c, that was obtained by Mekjian and Kapusta (15 



for kaons in a hadron gas at T = 180 MeV. Travelling outward in longitudi- 
nal and (more slowly) in transverse direction the average hadron encounters 
a rapidly diluting and cooling environment. In fact a combined analysis of 
pion pair correlation and hadron spectral data for central Pb + Pb SPS colli- 
sions [|14j] has revealed evidence for collective outward velocity fields reaching 



up to (3± ~ 0.55 (reminiscent of the limiting hydrodynamical fluid velocity 
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of c/y/3), the temperature falling to T = 120 MeV at hadronic decoupling. 
This fall thus occurs during an average time span of 6 fm/c, during which 
the total reaction volume increases about tenfold |Q [CJJ]. While the average 
hadron encounters about 6 binary collisions perhaps maintaining local 



momentum equilibrium, the inelastic fraction fades away rapidly in such an 



"explosive" hadronic expansion |T7|]. Thus the primordial population pattern 



created by the hadronization transition among hadronic and resonance species 
stays essentially unaltered (the resonances decaying within their proper time), 
remaining "frozen in" throughout the hadronic expansion phase. 

The fingerprint of the parton to hadron phase transformation in central 
Pb + Pb collisions is thus preserved in the observed bulk hadron yields, and 
yield ratios. Moreover the models of microscopic parton transport toward the 
critical QCD conditions of hadronization indicate statistical phase space dom- 
inance, according to the spin, mass and flavour spectrum of created hadronic 
and resonance states. Not surprisingly, the thermal hadrochemistry models 
[0, H [19| thus come up with meaningful accounts of this initial maximum en- 



tropy state. This lifts the conceptual qualms ever prevailing since Hagedorns 
engagement |24[] with statistical thermodynamics applied to hadron produc- 
tion data. There is no chemical equilibration occuring at the hadronic level 
to dynamically achieve a state of limiting temperature, neither in elementary 
e + + or p + p collisions nor in ultra- relativistic nuclear collisions |25j . Such 
states are created from partonic energy densities and temperatures above the 
limiting hadronic conditions, the QCD hadronization thus populating a sta- 
tistically ordered system (to avoid the term "state" as it invites the view of 
stationary conditions) just at the limiting hadronic energy density. This results 
in the apparent universality |7]] of hadron production as reflected in thermal 
model analysis from e + + e~ to Pb + Pb. 

In summary once we take for granted the tentative evidence from J/ty 
suppression, and from estimates of primordial energy density, in favour of the 
hypothesis that SPS Pb + Pb collisions reach beyond the limit of the hadronic 
phase of matter (as predicted by lattice QCD theory) the system has to re- 
hadronize upon expansion. And, second, the properties of this long-sought 
parton to hadron phase transition are reflected in characteristic features of 
the bulk hadronic phase that survive the rescattering and expansion phase 
because they freeze-in right at, or near the phase transition. Our focus was 
placed on the regularities of hadronic yield ratios because they are allowing 
for an estimate of the critical transition temperature (and the corresponding 
value of the baryo-chemical potential \ib that we did not comment about) 
owing to the applicability of thermal model analysis to this maximum entropy 
state that we have tried to establish. One now expects a convergence among 
the various thermal model approaches JT2], [8], IS], to finally pin down T and \xb 
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of the phase transition. At present we may conclude that the QCD transition 
temperature T c = 170 - 180 MeV, at \i B « 0.25 GeV. 

A final question requires attention: does universality of hadronization im- 
ply that small and large systems hadronize by the identical same QCD mech- 
anism? The author is not aware of a rigorous QCD argument. The colour and 
flavour composition of the partonic phase should reflect in the hadronization 
outcome (if there are mostly gluons, if strange quarks are far underpopulated 
or near-equilibrated etc.). Moreover there are rather straight forward effects 
of increasing system size. In the phenomenological coalescence approach El- 
lis and Geiger || have studied the step from e + + e~ — > Z° — > hadrons to 
e + + e~ — > W + W~ — > hadrons (at LEP 2). In the latter case, the two parton 
cascades may overlap. This merging of the decaying subvolumes leads to to- 
tal multiplicity increases, and the transverse momentum spectra exhibit some 
softening. The LEP 2 data have not been analyzed yet to demonstrate such 
effects. However, comparison of elementary p+p data with central Pb + Pb col- 
lisions shows a drastic effect in the hadronic strangeness composition: the ratio 
of strange to nonstrange quarks contained in the hadrons increases by a factor 



2 [12, 18], the ratio of double strange cascade hyperons to nucleons by a factor 



of about 10 [26] , the same ratio for ft, hyperons by a factor approaching 20 
[27|. Thus there is no flavour scaling from the elementary nucleon-nucleon to 
central Pb + Pb collisions. A certain fraction of this effect has been explained 
in the hadro-chemical model |12|, |28| to arise from the so-called "canonical 
suppression" (e. g. of strangeness) dominating the small system, due to local 
difficulty of energy and quantum number conservation [2£j. This effect disap- 
pears with increasing system size, a precursor seen in the overlapping W + W~ 
cascade decay mode ||, and plays no role in central Pb + Pb collisions which 
thus receive a grand-canonical treatment in these models. However, there is 
still further strangeness enhancement beyond this mechanism, which may be 
traced to before the phase transition assuming a higher relative strangeness 



content in the partonic Pb + Pb collision phase than in p + p collisions 
Primordial flavour or gluonic composition levels, perhaps characteristic of a 
symmetric partonic phase approaching the ideal QCD plasma state |28] may 



thus survive even the processes occuring during the hadronization phase tran- 



sition 1 30 1 . This then reflects in additional (e. g. flavour) parameters in the 
hadrochemical model. This new aspect (see, however, ref. ||) adds further 
significance to bulk hadronic signals from SPS nucleus-nucleus collisions. 
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Figure 1: Time profile of the pion decoupling strength from the hadronic fire- 
ball formed in central Pb+Pb collisions (arbitrarily normalized to unity). From 
two pion correlation analysis jTJ], |2T], we get an average decoupling time 
< Tf >= 8 fm/c, and a mean duration of the decoupling phase characterized 
by a width At = 4 fm/c for the assumed Gaussian profile. 



9 



1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 




Pb + Pb central v^s =17.2 GeV 






i y-y C M 1 < 05 






1 z 1 < 0.5 fm 






1 r_i_ 1 < 4 fm 






Geiger and Srivastava 




quarks and gluons 




clusters and hadrons 




II 1 1 1 1 1 1 1 


i 1 i i i i 1 i i i i 1 







10 
t (fm/c) 



15 



20 



Figure 2: Parton cascade model prediction for the time profile of the energy 
density contained in partons and hadrons, as found in a central cylinder sub- 
volume of 4 fm/c transverse radius at \y — Ucm\ < 0.5 in central Pb + Pb 
collisions at y/s = 17.2 GeV. Adapted from Geiger and Srivastava [f23|j . 
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